Quantum Criticality in Topological Insulators and Superconductors: Emergence of 
Strongly Coupled Majoranas and Supersymmetry 



(N 

o 

(N 
C 



i 

C 

o 
o 



> 

m 
m 

o 

(N 



X 



Tarun Grover 1 and Ashvin Vishwanath 1 
1 Department of Physics, University of California, Berkeley CA 94720 

We study symmetry breaking quantum phase transitions in topological insulators and supercon- 
ductors where the single electron gap remains open in the bulk. Specifically, we consider spontaneous 
breaking of the symmetry that protects the gapless boundary modes, so that in the ordered phase 
these modes are gapped. Here we determine the fate of the topological boundary modes right at the 
transition where they are coupled to the strongly fluctuating order parameter field. Using a combi- 
nation of exact solutions and renormalization group techniques, we find that the surface fermionic 
modes either decouple from the bulk fluctuations, or flow to a strongly coupled fixed point which 
remains gapless. In addition, we study transitions where the critical fluctuations are confined only 
to the surface and find that in several cases the critical point is naturally supersymmetric. This 
allows a determination of critical exponents and points to an underlying connection between band 
topology and supersymmetry. Finally, we study the fate of gapless Majorana modes localized on 
point and line defects in topological superconductors at bulk criticality, which is analogous to a 
quantum impurity problem. Again, an interplay of topology and strong correlations causes these 
modes to remain gapless but in a strongly coupled state. Experimental candidates for realizing these 
phenomena are discussed. 

PACS numbers: 



The experimental verification of Z2 topological insula- 
tors in materials with strong spin-orbit coupling [H-0] has 
lead to an explosion of interest in topological phenom- 
ena in solids. The theoretical description of these phases 
is now well understood, since the key properties such as 
the protected edge states are captured by models of non- 
interacting fermionsQ- 17 1. Experimental realizations of 
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Figure 1: The two class of problems studied in this paper 
that concern the fate of fermionic surface modes in TI/TSC at 
a critical point. T denotes temperature while r is the tuning 
parameter for the quantum phase transition. The counter- 
propagating fermionic surface modes (blue lines) interact with 
the critical fluctuations (headed vectors). The critical point 
may correspond to a bulk transition as depicted in (a) or it 
could be a surface transition where the bulk is unaffected as 
shown in (b). 



these states in various bismuth based materials typically 
feature weak electron correlations, and are well described 
by conventional band theory. However, recent work has 
begun to study the interplay of electron correlations and 
band topology, either by introducing correlated transi- 
tion metal ions such as Mn, Fe and Cu in bismuth based 
topological insulators [Tel . 0] , or by studying correlated 
materials with strong spin-orbit interactions such as irid- 
ium based oxides and f-electron systems US-El. Numer- 
ical studies of model correlated systems such as Kane- 
Mele-Hubbard Hamiltonian 0, HH have also appeared. 
An important question for theory then is - are there qual- 
itatively new phenomena that emerge from strong corre- 
lations, that are not captured by perturbing free electron 
topological states? 

A novel feature of Z2 topological insulators, in con- 
trast to integer quantum Hall states, is that they are well 
defined only in the presence of time reversal symmetry. 
A complete classification fl4l . [lit °f f ree fermion topo- 
logical phases reveals a superconducting analog of these 
states, which akin to topological insulators lack spin ro- 
tation symmetry, and which are also only well defined 
only in the presence of time reversal symmetry (class 
Dili). The well known B phase of superfluid He-3 is a 
realization of a topological phase in this class (27I I32I - 
[35|, and possible solid state realizations have also been 
proposed 28, 36, 3^1 . In such symmetry protected topo- 
logical phases, spontaneous symmetry breaking provides 
a new route to exit the topological phase, without clos- 
ing the fermionic energy gap in the bulk. For example, 
electron correlations often lead to magnetic order, which 
breaks time reversal symmetry. A natural question con- 
cerns the evolution of the topological phase as magnetic 
order sets in. Therefore we consider quantum phase tran- 
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sitions from a paramagnetic topological phase, protected 
by time reversal symmetry, to a magnetically ordered 
phase, where the topological surface states acquire a gap. 
Although critical properties of the bulk are unaffected by 
the band topology, critical fluctuations near the surface 
will interact with topological surface modes. Note that 
these magnetic fluctuations have a significant impact on 
the protected surface modes when they condense, lead- 
ing to a gap. Do the surface states survive at the critical 
point, and if so are they fundamentally different from 
the free fermion states of the paramagnetic phase? Since 
critical fluctuations are strongly interacting, particularly 
in low dimensions, answering this question typically re- 
quires non-perturbative techniques. We use a combina- 
tion of exact solutions, renormalization group and map- 
ping to models with enlarged symmetry, to study class 
Dili topological superconductors (DIII-TSC) in d=l,2,3 
spatial dimensions and topological insulators in d=2,3. 
As usual, our results apply to a finite temperature 'quan- 
tum critical' fan of the phase diagram, which is controlled 
by the zero temperature critical point (38|. 

Our main results are: (i) At the bulk critical point, 
the surface modes of the d=l time reversal symmetric 
(Dili) TSC, which are a pair of Majorana zero modes, 
remain gapless but couple strongly to the critical fluc- 
tuations. This leads to edge fermion correlation func- 
tions that depart qualitatively from that of free Majo- 
rana fermions. We calculate their precise form using an 
exact solution, (ii) However, for DIH-TSCs in d=2,3 or 
TIs in d=3, the gapless surface modes decouple from the 
critical fluctuations at low energy, and are qualitatively 
unaffected at the critical point. These results depend cru- 
cially on the fact that surface exponents of a transition 
are rather different from those in the bulk, (iii) We also 
consider the case when the surface orders magnetically 
before the bulk. In particular we study the boundaries 
of DIII-TSC in d=2,3, when magnetic order sets in on 
crossing a quantum Ising critical point. These systems 
are found to display the remarkable feature of emergent 
supersymmetry at the critical point. The gapless criti- 
cal modes of the Ising transition are superpartners of the 
Majorana surface states. In d=2 we argue that this is the 
tricritical Ising transition, accessed however by tuning a 
single parameter to attain surface criticality. Similarly, 
the d=3 system is argued to realize a supersymmetric 
critical point by tuning just a single parameter to at- 
tain surface criticality. Since this is one of the simplest 
settings to realize supersymmetry in a condensed mat- 
ter system, we also discuss a number of striking physical 
consequences. A similar phenomena occurs at the surface 
of TIs in d=3, at the critical point into a paired super- 
conductor, which however requires tuning more parame- 
ters (doubled versions of this were considered previously 
in different contexts (39l . Finally (iv) we discuss 

the fate of topologically protected modes at the cores of 
point and line defects, at the critical point. These open 



up a class of problems analogous to quantum impurity 
(or Kondo) problems, coupled to bosonic/fermionic bath 
[4l|-|46( , where the defect modes play the role of the impu- 
rity, coupled to bulk critical fluctuations. However, these 
defects could also be one dimensional modes propagat- 
ing along the defect line and therefore, apart from be- 
ing of interest to the physics of TI/TSC, these problems 
also serve to generalize the quantum impurity problem 
to higher dimensional 'impurities', that are now one di- 
mensional objects. For DIII-TSC we find that the point 
(in d=2) and line defects (in d=3) remain gapless at the 
critical point but are essentially modified from the free 
Majorana form. The precise form of the strong coupling 
theory in d=2 is accessed using an e(= 3 — d) expansion. 

We briefly review related prior work. Quantum phase 
transitions from a TI to a metal or trivial insulator, while 



preserving symmetries were studied in Ref. [47l . l48| . These 
transitions are rather different from the one described 
here, since they involve an intermediate state with gap- 
less delocalized fermions in the bulk. In the problems 
we consider, the fermions always remain gapped while a 
bosonic symmetry breaking order parameter condenses. 
Transitions where the symmetry protecting the topologi- 
cal phase is spontaneously broken have also been consid- 
ered previously, but their effect on the surface modes has 
not been discussed. Purely surface magnetic transitions 
in TIs were discussed by Xu [49$ . but the surface mag- 
netic phase transitions in DIII-TSC, which turn out to 
have emergent supersymmetry as we show in this paper, 
were not considered. Surface superconducting transitions 
in 3dTIs, which also turn out to be supersymmetric when 
chemical potential is also tuned to obtain Dirac point, 
have also not been discussed earlier. Quantum critical- 
ity in the Kane-Mele Hubbard model has been discussed 
in Refs.dSHl but again the impact on surface modes 
have so far been ignored. We now begin with a summary 
of the models and other preliminaries required for the 
remainder of the paper. 



MODELS AND MAIN RESULTS 

In this section, we provide a brief summary of the sys- 
tems studied in this paper along with a synopsis of our 
main results. Topological insulators (TI) and topologi- 
cal superconductors (TSC) are recently discovered sys- 
tems where the topology of underlying electronic band- 
structure leads to protected surface modes @. There 
are five different classes of TI/TSC in every dimension 
where different classes differ from each other based on the 
underl yin g symmetry (time reversal, particle-hole etc.) 
[II fl6l.Tl7|. Each class can itself contain a countably in- 
finite number of distinct phases ('Z classification'), two 
phases ('Z2 classification') or only one distinct phase ('Zi 
classification'). We focus on two classes: All topologi- 
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cal insulators (or simply topological insulators) and Dili 
topological superconductors, both of which require time- 
reversal symmetry. All has the classification Zi , Z 2 , Z 2 in 
dimensions d = 1, 2, 3 respectively while Dili is classified 
as Z2, Z2, Z in d = 1,2,3 respectively. 

We are primarily interested in the fate of bound- 
ary modes in these particular TI and TSCs at a mag- 
netic and/or superconducting symmetry breaking quan- 
tum phase transition, where the gap to the single-electron 
excitations in the bulk remains non-zero at the transition. 
Therefore, let us briefly discuss the nature of boundary 
states in these topological phases away from the phase 
transition. 

First consider All TIs. In d = 1, the classification 
is Zi and hence there are no protected edge modes. In 
d = 2, the canonical example of an All TI is a quantum 
spin-Hall insulator and a generic edge is described by 
an interacting helical Luttinger liquid [El], 55 1. In d = 
3, the All class corresponds to Z 2 topological insulators 
whose surface modes are given, in their simplest for m, by 
a single species of non- interacting Dirac fermion llMl3] . 



Weak interactions are irrelevant for these surface modes. 

In contrast, the Dili TSC have non-trivial edge modes 
in all three dimensions. In d = 1, these edge modes 
correspond to a pair of Majorana modes localized at the 
two ends of the system, that are related to one another by 
time reversal symmetry [34]. In d = 2, they correspond 
to a non-interacting helical Majorana liquid. Finally, in 
d = 3, again one obtains two-dimensional non- interacting 
Majorana modes. The integer classification in d = 3 is 
related to the minimum number of such Majorana surface 
modes required by bulk topology. Here we only consider 
the case with a single surface mode. In all three cases, 
the edge modes remain gapless and non-interacting for 
weak interactions. 

Strong interactions can lead to two very different kinds 
of symmetry breaking phase transitions in these topo- 
logical phases. In the first kind, the transition occurs 
throughout the system. In the second kind only the sur- 
face undergoes phase transition while the bulk remains 
gapped. We find that these two kinds of transitions com- 
paratively have very different influence on the fermionic 
edge modes. In the former case, the nature of critical 
fluctuations near the boundary of the system is dictated 
by the so-called 'boundary critical phenomena' (56, 57|. 
The primary implication of boundary critical phenomena 
is that the universal critical exponents associated with 
various quantities such as specific heat are rather differ- 
ent for the boundary fluctuations when compared with 
the bulk. This fact has important consequences for the 
fermionic edge/surface modes in TI/TSC at criticality. 
We find that in both two and three dimensions, owing to 
the non-trivial boundary critical exponents, the fermionic 
modes at the boundary completely decouple from the 
critical fluctuations for weak coupling. Therefore one ob- 
tains non-interacting gapless fermionic boundary modes 



even at the criticality. Only in the case of d = 1 Dili class 
TSC does the bulk fluctuations modify the behavior of 
fermionic edge modes. In this case, we find that even 
though the Majorana edge modes Xi>X2 couple strongly 
with the critical fluctuations, they still remain gapless 
and their density of states N(E) changes from being a 
delta function N(E) ~ 5(E) away from the criticality, to 
N(E) ~ constant at the critical point. Since the elec- 
tron creation operator is related to Majoranas through 
S = xi + «X2, the knowledge of N(E) allows one to cal- 
culate the local density of states of the electrons N e (E) 
at the end of the chain, which is precisely what would be 
measured by scanning tunneling microscopy (STM): 



N e (E) = Im / dt (c(r = 0,t = 0)c f (r = 0, t))e^ E+i0 ^ 



The differential conductance dl/dV in STM is given by 
dl/dV (x N e (E — eV) where V is the voltage difference 
between the tip and the sample and / is the current. 

The critical exponents associated with magnetic corre- 
lations C(r,t) = (S(r, t) S(0, 0)) near the boundary can 
in principle be accessed by Nuclear Magnetic Resonance 
(NMR) from the relaxation rate 1/Ti, which measures 
the imaginary part of spin-susceptibility \ a t tempera- 
ture T: 



^ = Limit (w -> 0) J d 



d k Imx(fc, uj) 



Imx is proportional to the fourier transform of the 
correlation function C (r, t) . 

For the second kind of phase transitions, where only 
the surface undergoes phase transition, the fermionic 
modes live in the same dimension as the critical fluc- 
tuations and the 'boundary critical phenomena' is not 
relevant. Correspondingly, we find that the boundary 
fermionic modes are strongly coupled to the critical fluc- 
tuations at the phase transition in almost all cases. Per- 
haps most interestingly, we find that an Ising magnetic 
phase transition in a Dili TSC has an emergent super- 
symmetry at criticality at long wavelengths without any 
additional fine-tuning in both two and three dimensions! 
We recall that the conventional symmetries of a quan- 
tum mechanical system, such as a rotation or a transla- 
tion, transform a boson into boson and a fermion into 
fermion, the generator of the symmetry thereby being a 
boson. The hallmark of supersymmetry is the invariance 
under rotation matrices that are themselves fermionic in 
nature. Therefore, supersymmetry rotates a fermionic 
operator into a bosonic one and vice- versa. This has im- 
portant physical consequences for our system, namely, 
Dili class TSC at criticality. For example, it results in 
the equality of anomalous dimension corresponding to 
the magnetic order parameter with that of the Majorana 
fermion, precisely because the bosonic order parameter 



4 




HtUH 



h-h. 



(0,0) 



Figure 2: (above) Illustration of Majorana boundary modes 
771 ,772 in a one dimensional time-reversal invariant TSC cou- 
pled to critical Ising fluctuations with interaction strength g. 
(below) Renormalization group flow for the same problem, h 
is the transverse field for the Ising magnet. The g = fixed 
point is unstable to g* = 1 fixed point where Majorana modes 
are strongly coupled to the critical fluctuations and have cor- 
relation function (g a (0)g a (r)) ~ 1/r for a = 1, 2. 



can be rotated into a fermionic Majorana via supersym- 
metry. We also find emergent supersymmetry at a su- 
perconducting phase transition on the surface of a three- 
dimensional All TI with chemical potential tuned so as 
to obtain a Dirac fermion coupled to critical supercon- 
ducting fluctuations. 

Dili TSCs also host gapless Majorana modes at certain 
point and line defects in d = 2 and d = 3 respectively 



|3J, l58| . The Hamiltonian for these modes in the absence 
of any interactions is identical to that of the edge modes 
in d = 1 and d = 2 Dili TSC. We study the fate of 
these Majorana modes at bulk criticality and find that 
the point defect exhibits anomalous exponents, and is 
hence very different from a free Majorana mode. The 
line defect behaves like a 'marginal' Majorana liquid at 
criticality, with logarithmic corrections to free Majorana 
correlation functions. 



SURFACE FERMIONIC MODES COUPLED TO 
BULK CRITICALITY 

(i) Majorana modes in one-dimensional TSC: 
We begin by considering the simplest example in this 
class of problems: a one-dimensional TSC in class Dili 
at the criticality between TSC and an Ising magneti- 
cally ordered ordinary superconductor. For simplicity we 
consider a semi-infinite wire of a TSC and consider the 
case of a finite chain in the Supplementary Information, 
Sec. Ai The Hamiltonian of the system may be written 
as: 



H — Hjsing -\- H I sing — Majorana 
00 



where xi and Xi are the two Majorana modes located 
at the end of the wire. The fermionic degrees of the free- 
dom are gapped in the bulk on either side of the phase 
transition and do not play any role in the following dis- 
cussion. Therefore we have not included them in the 
Hamiltonian H (this holds for all phenomena considered 
in this paper). When h > 1, the Ising order parame- 
ter (a z ) = and the coupling g between the Majorana 
modes and the Ising field erf is irrelevant, that is, there 
are (free) Majorana modes at the end of the wire, as 
they should, in a TSC. On the other hand, when h < 1, 
(a z ) 7^ 0, which provides a mass to the edge Majorana 
modes. We are interested in the fate of the Majorana 
modes right at the criticality (h = 1). Using a clas- 
sic result from the theory of boundary conformal field 
theory (BCFT) [Hoi . [6l| , the scaling dimension for the 
boundary order spin a\ is 1/2. Therefore the interaction 
term gi J dt Xi{ T )X2ij)a z (x = 0, r) between the Majo- 
rana spins and the spin-chain is relevant at the decou- 
pled fixed point (g = 0) in the renormalization group 
(RG) sense and hence it is an unstable fixed point. To 
proceed, let us define Cq = 1x1X2- Further, let us set 
g = 1. We will soon consider the deviation of g from 
unity. Thus, at criticality, H in Eqn|T] may be rewritten 



H = - 



OO 

E 

=1 



i+1 



(2) 



i=0 



H as written above corresponds to transverse field 
Ising model without a transverse field for the first spin 
(To- Following [iH], let us perform the following duality 
transformation: 



Mo 

M f_iM| V i > 1 
/4 V i > 



(T4 & 



1+1 



Hamiltonian H in terms of the new variables fi is given 
by the standard transverse field Ising model: 



H 



E 



Mi + N n\ 



+1) 



This form allows one to use the results from boundary 
conformal field theory to calculate the unequal time cor- 
relation function for Majorana fermions xii X%- I n partic- 
ular, the spin-operator (Tq may be denned as ctq = iX3Xi 
where X3 is a Majorana operator without any dynamics. 
This definition ensures that Cg anticommutes with ctq. 
Therefore, the time-correlation functions of %i are same 
as that of operator <7g (= /jg): 



(Xi(O)xi(t)) = </^(0K(t)) ~ 1/r 
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Figure 3: The coupling to bulk critical fluctuations is ir- 
relevant for two and three dimensional Dili TSCs and three 
dimensional All TIs. This means that at low energies, the 
boundary fermionic modes completely decouple from the bulk 
critical fluctuations. 



since the scaling dimension of the boundary spin 
is 1/2 in Ising BCFT. Due to the symmetry, the cor- 
relation functions of \2 are identical to that of %i- 
Let us consider deviation of g from unity: the term 
(g — l)«XiX2 (J o = (s — 1)^0 turns out to be irrelevant 
because the boundary scaling dimension for the trans- 
verse field is 2 [6lJ. Therefore, g = 1 is a stable fixed 
point. We conclude that the interaction with the critical 
modes dramatically changes the scaling dimension of the 
Majorana modes from to 1/2 and therefore they behave 
like a zero-dimensional non-Fermi liquid. The density of 
states at energy E changes from being a delta function 
5(E) in the non- interacting limit to a non-zero constant 
at the stable fixed point which can be probed using STM. 
The feedback of the Majorana fermions on the Ising chain 
is also rather dramatic. When g = 0, the edge correla- 
tions (of (0)crf (r)) ~ 1/t, as discussed above. Using the 
above duality, at g = 1, that is, at the stable fixed point, 
the operator a\ = <Jq(J.q. Since, ctq commutes with the 
Hamiltonian and therefore does not have any dynamics, 

(of (°)°i( r )) = (Mg(0K( T )) ~ V^ 4 - This follows from 
the fact that the scaling dimension of the energy oper- 
ator at the boundary (= /ig) is 2 j(3l||. This exponent 
can be measured from the imaginary part of the spin- 
susceptibility. The renormalization group flow diagram 
is shown in Fig|2] 

(ii)Majorana modes in two/three-dimensional TSC: 
Let us study analogous problems in higher dimensions. 
As a first example, consider the magnetic phase transi- 
tion out of a two dimensional TSC in class DHL The un- 
perturbed edge states consist of helical Majorana modes: 



H = J dx iv F (xrOxXr - XlOxXl) 



(3) 



At the critical point, these edge states are coupled to 
the critical magnetic fluctuations <j) via a term AS in the 
action 



AS = ig / dx\\dTXL(T,x\\)XR{T,x\\)(l>(T,x±. =0,Z||) 



where xm(x_l) denote spatial coordinates parallel (per- 
pendicular) to the boundary of the system. To deter- 
mine the relevance/irrelevance of this term at the g = 
fixed point, similar to the case of Majorana mode above, 
it is crucial to use the boundary critical exponents for 
the magnetic fluctuations 0. The relevant exponents are 
known to a very good accuracy numerically [62j . In par- 
ticular, the edge correlation functions (4>(x± = 0, xii = 
0)<f>(x± = 0,X|| = L)) decay approximately as [g^ 



((j)(x ± = 0,x\\ = 0)(/){x ± = 0,x,| = L)) ~ 1/L 



2.r>2 



(5) 



The non-interacting Majorana fermions in two space- 
time dimensions have a scaling dimension of 1/2. There- 
fore, in renormalization group sense, the coupling g at 
length scale L scales as g(L) ~ 1/L 26 making it ir- 
relevant at the 5 = fixed point. Therefore, for weak g, 
one obtains well defined non-interacting Majorana modes 
even at the critical point between the TSC and a mag- 
netically ordered trivial superconductor. Note that if one 
naively uses the bulk scaling dimension for 4>, which is 
close to 1/2, one would instead reach the erroneous con- 
clusion that the coupling g is relevant. This again under- 
scores the importance of boundary critical phenomena for 
the fate of topological surface modes at the bulk ordering 
transition. 

Next, consider the same problem for a three dimen- 
sional Dili TSC. The boundary fermionic modes in this 
case are described by the Hamiltonian 



H 



= I d 2 x iv F X T (<?zd x 



xdy) X 



(6) 



(4) 



where x = [Xi X2] is a two-component Majorana 
spinor. The coupling between these surface modes and 
the critical fluctuations is given by an expression similar 
to Eq0 



AS = ig J d X\\dTXL(T,X\\)xR(T,X\y)(l>(T,X_L=0,X\\) 

(7) 

where xu is now two-dimensional. In this case, the 
space-time dimension for the bulk is four, which is the 
upper critical dimension for the critical fluctuations and 
therefore the boundary critical correlations exhibit mean- 
field behavior |H(||57j: (4>( x ± = 0, aril = 0)4>(x±_ = 0, xu = 
L)) ~ 1/L 4 . The scaling dimension of Majorana fermions 
in three space-time dimensions being unity, this implies 
that g is again irrelevant. Thus we find that the surface 
Majorana modes remain well-defined for small values of 
coupling to the critical modes for a three dimensional 
superconductor as well. 

(iii) Dirac modes in three dimensional TI: 
The analysis of three dimensional All TIs where 
the chemical potential has been tuned to obtain Dirac 
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fermions at its surface proceeds in the exact same man- 
ner as that for three dimensional Dili TSCs and one 
finds that Dirac fermions remain well-defined at the crit- 
icality. Note that if the critical fluctuations were instead 
confined to the boundary (that is, a surface phase tran- 
sition), the result would be entirely different. In particu- 
lar, for a surface transition the coupling to critical modes 
will be relevant and the Dirac fermions will not be free 
at the transition [ijj- For a generic chemical potential 
one obtains a helical Fermi liquid. Ref.(49| studied the 
problem of helical Fermi surface at a surface phase tran- 
sition and it was found that the coupling between helical 
Fermi liquid and critical fluctuations is irrelevant. For 
bulk critical fluctuations, the problem of interest in this 
section, the Landau damping term has the same form 
as in Ref.[49| ~ J dujd 2 q\\dq±\4>(q\\, q±, uj)\ 2 \uj\q\\ and the 
coupling to critical fluctuations is even more strongly ir- 
relevant owing to the boundary critical exponents cited 
above. Therefore, the helical Fermi liquid is unaffected by 
bulk critical fluctuations for weak coupling. The case of 
two dimensional TIs and the Kane-Mele Hubbard model 
requires a more careful analysis due to the Luttinger liq- 
uid nature of the edge states and the results will be re- 
ported elsewhere. 



SURFACE FERMIONIC MODES COUPLED TO 
SURFACE CRITIC ALITY 



In the previous section we studied the fate of fermionic 
boundary modes at bulk symmetry breaking quantum 
phase transitions. However, it is also conceivable that in 
certain situations the surface undergoes transition before 
the bulk does. In this section, we study the fate of surface 
modes at such surface phase transitions. 

(i) Emergent supersymmetry in three-dimensional TSC: 
Let us consider a surface magnetic phase transition in 
a three dimensional Dili TSc. The unperturbed surface 
states are described by the Hamiltonian in Eqn. [Bl At 
the transition, these surface states are coupled to an Ising 
magnetic order parameter. Note that one cannot couple 
the system to an XY or 0(3) magnetic order parameter 
since the discrete time-reversal is the only symmetry that 
the surface states possess. A scenario such as this can 
arise in several ways. One can imagine a situation where 
the material has two bands, one of them contributes itin- 
erant electrons which form a topological superconductor 
while the other band consists of localized moments that 
provide the Ising magnetic variable. Alternatively, one 
can also dope the system with magnetic impurities in a 
uniform manner so that the disorder effects can be ig- 
nored. What is the fate of the boundary states at such a 
surface phase transition if the bulk still remains gapped? 
The low-energy theory of the Majorana modes coupled 
to Ising critical fluctuations is: 



S = J drdxdy [\ T (a y d T + a x d x + a z d y ) \ + ig 4> X T X+ 
I 



(8) 



where x = [xi X2} io'y If 9 were zero, the Ising or- 
der parameter will undergo a phase transition such that 
the critical point is described by the conventional Wilson- 
Fisher fixed point [63j |. However, as shown in Ref . |64 . 65 1. 
in the presence of the coupling g, the Wilson-Fisher fixed 
point is unstable to a new conformally invariant fixed 
point that has an emergent M = I supersymmetry. As 
mentioned earlier that supersymmetry allows one to ro- 
tate fermions into bosons (and vice-versa) and the label 
Af denotes the independent number of such fermionic ro- 
tation operators. Supersymmetry requires u = g 2 /8 and 
one would have naively thought that satisfying this re- 
lation requires fine tuning and hence the supersymmet- 
ric fixed point would be multi-critical. However, in this 
case, the supersymmetry is emergent (FmI. ItiBl). by which 
we mean that this relation is satisfied automatically at 
the fixed point and the only tuning parameter is the bo- 
son mass r, as in a conventional phase transition. Let 
us briefly consider the supersymmetric transformations 
that keep the above action invariant (o^ : 
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$X 
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ei £2] T is 



where 7^ = {cr y ,a x ,a z } and e 
a two-component Majorana variable. The emer- 
gence of this 'fermionic rotational invariance' implies 
that at the critical point the correlation function of 
the Majorana fermions (x T (w,p)x(— w, — p)) 



p 



{(JyUi + a x p x + o z py)~ l and the order parameter 
{(j>(oj,p)4>(—uj,—p) ~ (p 2 + lo 2 ) 71 / 2 ^ 1 have exactly the 
same scaling exponent rj (65l |. Experimentally, this can 
be verified by extracting the fermionic anomalous expo- 
nent using STM or ARPES while that corresponding to 
the order parameter can be extracted from the neutron 
scattering. 

At finite temperature, the supersymmetry is sponta- 
neously broken (o7l - l69j | and as a result, one obtains a 
collective fermionic excitation, 'phonino', an analog of 
Goldstone boson, as first pointed out by Lebedev and 
Smilga [7(i| . Since the physics in quantum critical regime 
is controlled by the critical fixed point, the contribu- 
tion of phonino to various thermodynamic and transport 
properties should be visible in this regime. We leave the 
details of this interesting subject to future work. 

(ii) Emergent supersymmetry in two-dimensional TSC: 
Let us consider a one lower dimensional analog of the 
problem studied above viz. boundary of a two dimen- 
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sional Dili TSc (Eqn(3]) coupled to critical Ising magnetic 
fluctuations: 



S = J drdx [{xiidrXR + XhdrXL + iXudxXR - iXhdxXL 
{\{d^? + \<? + j^) + KgxLXR^ 

In the absence of any coupling to the magnetic fluctua- 
tions, the edge theory of the Majorana modes is identical 
to that of a 1 + 1-d critical Ising model self-tuned to crit- 
icality. Therefore when 5 = and the <f> field is at Ising 
criticality, the above action describes two decoupled crit- 
ical Ising models. The scaling dimension A g of g equals 
A g = 1 + 1/4 = 5/4 < 2, and therefore it is a relevant 
perturbation at this decoupled fixed point. What is the 
fate of the coupled theory? Below we present several ar- 
guments that it flows to the tricritical Ising fixed point 
with conformal charge c = 7/10. 

Let us begin with a simple mean field analysis and 
denote the Ising variable corresponding to the bosonized 
Majorana modes as (j>\ and define <p 2 = </>• At the mean 
field level, the coupled system may be described by the 
following Landau free energy: 

F = ui4>\ + vx4>l + r 2 4>l + u 2 4>i + v 2 <t>t + g<t>l4>2 (10) 

<f>i term in F is not allowed since the Majorana modes 
are self-tuned to criticality, being the boundary of a TSC. 
Furthermore, note that the g term breaks the symmetry 
from Z 2 x Z 2 to Z 2 as in Eqn(5] and describes an energy- 
order parameter coupling. Assuming <f> 2 to be small near 
the critical point, one can eliminate it in favor of (f>\ 

,2 

<fi 2 « — J^r- • This leads to the following free energy (upto 
0(4>\)) for the variable <f>\: 

F = (u 1 -^-)4>i + v 1 cf> 6 1 (11) 
2r 2 

When r2 > 1, the minimum of F lies at <pi — with 
F vanishing as quartic power of (f>\. This implies that 
the Majorana fermions are at Ising criticality while the 
field 4> is gapped consistent with the action in Eqj9j At 
r 2 = l/2u\, the potential for </>i becomes ~ <f>\ which cor- 
responds to Ising model at tricriticality. For r 2 < l/2u\, 
the minima of F lie at 4>\ + 1 which corresponds to mass 
acquisition of Majorana due to condensation of <fi. There- 
fore, we conclude that within the mean-field theory, the 
phase transition from a non-magnetic one- dimensional 
TSC to a magnetic superconductor is described by a tri- 
critical Ising model! 

There are several consistency checks that support this 
conclusion. First, Zamolodchikov's c-theorem |7l[ im- 
plies that the central charge of the system in EqJS] at 
the magnetic phase transition must be less than that 



of two decoupled Ising models (= 1/2 + 1/2 = 1). In- 
deed, the central charge for the tricritical Ising model 
is c = 7/10, consistent with this constraint. Further- 
more, if the magnetic phase transition is indeed in the 
tpicrital Ising universality, owing to the c-theorem, a rel- 
evant perturbation can either take it to a gapped fixed 
^int (c = 0) or Ising critical point (c = 1/2), since those 
are the only two allowed values of central charge for uni- 
tary theories that are less than 7/10. This property is 
indeed shared by the action in Eq(9] giving a mass to <f) 
takes the system to Ising critical point while condensing 
<j) takes it to a gapped phase. Interestingly, the tricrit- 
ical Ising has emergent M = 1 supersymmetry (72I |73| , 
consistent with the matter content of the action in Eqn|!|] 
which contains one real fermion and one real boson. The 
Majorana fermion x an d the scalar <fi 2 are superpartners 
of each other with the same anomalous dimension of 1/5 
(6l) . Like in other cases, the experimental predictions 
such as the scaling dimension of Majorana fermion and 
the Ising order parameter can in principle be accessed via 
single particle tunneling and NMR. 

Interestingly, the flow from the supersymmetric Tri- 
critical Ising model to the Ising model can be thought of 
as spontaneous breaking of the supersymmetry. Within 
this picture, the free Majorana fermion at the Ising crit- 
ical point can be thought of as the 'goldstino' associated 
with the spontaneous breaking of the supersymmetry. 
This picture also generalizes to the higher dimensional 
supersymmetric models that we study. Since the pres- 
ence of free fermionic boundary modes in the topologi- 
cal phase is a consequence of topological band structure, 
which is well-defined even at the surface phase transition, 
this observation hints at an interesting interplay between 
band topology and emergence of supersymmetry at the 
surface quantum critical points. 

(iii) Emergent supersymmetry in three-dimensional TI: 
Next, we study surface phase transitions in a three 
dimensional All TI. In the absence of interactions, the 
Hamiltonian for a surface with normal ft = z is given by: 

H = ^2 ct (k x <r y - k y o- x - n)c (12) 

where /1 is chemical potential. Let us consider fine- 
tuning /i to zero so as to obtain fermions with massless 
Dirac dispersion at the surface and study the instabil- 
ity of these surface modes to an s-wave superconductor. 
Such a transition can potentially be driven by intrinsic 
interactions and may also be realized by pattering the 
surface with a Josephson junction array. We further re- 
strict ourselves to a case where the superconducting fluc- 
tuations are particle-hole symmetric and therefore, in the 
absence of coupling to the surface fermions, would ex- 
hibit a quantum phase transition with dynamic critical 
exponent z = 1 (this may require additional fine-tuning) . 
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Thereby the effective action for the coupled system near 
the phase transition may be written as 



(13) 



S = J d 3 x [tp^d^ip + g ip T eip + c.c.) 
+ |^| 2 +rH 2 + U |0| 4 l 



where (f> is the superconducting order parameter, 
i' T = [cf cj.], {70,71,72} = {<Tz,Vx,Vy} and e is two- 
dimensional antisymmetric tensor. Akin to the case of 
Dili TSC coupled to Ising order parameter, the above ac- 
tion is also known to flow to a supersymmetric fixed point 
[iol . 64 1 where the fixed point value of the couplings sat- 
t. In this case, the supersymmetry corresponds 



isfy g 2 -- 

to a N — 2 Wess-Zumino model which implies that there 
exist two fermionic generators that implement supersym- 
metry. Supersymmetry allows one to calculate the exact 
anomalous dimensions at this fixed point. In particular, 
7/0 and ?70 both equal 1/3 [75j . To further explore the 
experimentally testable consequences of supersymmetry, 
we note that like any other symmetry, supersymmetry 
is associated with conserved currents ('supercurrents'), 
which in this case are given by (66|: 



J 



- {lpd P 4>*) + i 2 < ^* ^ (*72)7m 



Being conserved, supercharges Jo, Jo do not ac- 
quire any anomalous dimensions at the critical point 
and therefore their correlation functions are given by: 
(J (t,x)J (t,x)) ~ ( X 2+ T i)2 ■ Tlie supercharge J car- 
ries electromagnetic charge of three units and trans- 
forms in the same way as the microscopic operator O ~ 
Of(x)c^(x)cf(x') where the points x and x' are close to 
each other. It might be possible to measure the correla- 
tion functions of O in an ARPES experiment similar to 
that used to measure cooper-pair correlations (76| . 

For completeness, let us briefly mention the results for 
magnetic phase transitions at the surface of a three di- 
mensional All TI, as dicussed in Ref . (49| . When the 
chemical potential is tuned to obtain a Dirac fermion, 
though the theory is not supersymmetric (recall that 
supersymmetry requires equal number of bosonic and 
fermionic species), the fermions are still strongly cou- 
pled to the Ising order parameter and acquire a non- 
trivial anomalous dimension, calculable in an e-expansion 



49. 



78j. For a generic chemical potential, one obtains a 
helical Fermi surface for the Dirac fermions and for weak 
coupling, the interaction between the fermions and the 
critical modes is irrelevant [4j| and the magnetic transi- 
tion remains in the three-dimensional Ising universality 
class. 



(iv) Surface transitions in two dimensional TI: 
In contrast to the above examples, the surface phase 
transitions in a two dimensional TI coupled to a power- 
law XY magnetic or superconducting order are of qual- 
itatively different nature. This is because the one- 
dimensional edge of a TI can support a stable power- 
law Berezinskii-Kosterlitz-Thouless (BKT) phase [79( at 
zero temperature and therefore, one can potentially find 
a whole phase (rather than a fixed point) where the heli- 
cal fermions and critical bosonic fluctuations are strongly 
coupled. As a first example, consider the case of helical 
edge modes of TI 55| coupled to critical magnetic fluc- 
tuations. The bosonized Hamiltonian is: 



H = J (kn uK x (d x 9) 2 + -jL {d^f + K 2 (d x £) 2 + 

-}-{d x v) 2 +r v cos{v) +0cos(2«f>-£) (14) 

Here we have bosonized the helical Luttinger liquid 
with ip m = e 1 ^-^, Vli = e i{ - e+ ® and [<p(x), V0(V)] = 
iirS(x— x'). u is the velocity of helical modes, £ is the XY 
magnetic order parameter field while e w inserts a 2ir vor- 
tex in the £ field and satisfies the commutation relation 
[£(x), Vf(x')] = iir5(x — x'). Depending on the magni- 
tude of K2, the vortex term r v can be relevant /irrelevant. 
When r v is relevant, the power- law XY-order is lost and 
the correlations of the (j> ne ld decay exponentially. In 
this case, the coupling g to the helical modes is irrele- 
vant and thus the fermionic edge modes are not affected 
by the critical fluctuations. On the other hand, when r v 
is irrelevant, the field x is m the BKT phase and one can 
always find K\,K 2 such that g is relevant. Relevance of 
g implies that the fields 2cf> and £ lock together and are 
described by a single compact variable. Therefore the 
central charge of the system jumps from two to one as 
g becomes relevant. The analysis for helical Luttinger 
liquid coupled to superconducting fluctuations proceeds 
analogously with the replacement <fi — > 9 in the interac- 
tion term g. The superconducting phase locks to the 28 
and the central charge jumps across the phase transition 
from two to one, gapping out half the degrees of freedom. 



LINE/POINT DEFECTS COUPLED TO BULK 
CRITIC ALITY 



Lastly, we discuss the fate of line and point defects 
in d=2 and d=3 Dili TSCs respectively, when the sys- 
tem is driven through a magnetic phase transition. Since 
the defects should preserve time reversal symmetry these 
cannot be regular vortices. Instead, a time reversal in- 
variant vortex which winds in spin space is required as 
discussed in [34[. Alternately, a lattice dislocation in a 
system with a weak topological index leads to gapless 
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Total spatial 
dimension d 4- 


Surface Modes at 
Bulk Criticality 


Topological Defects at 
Bulk Criticality 


Surface-only Criticality 


d = 1 


TSC Edge modes 
(x(r)x(O)) ~ 1/r 






d = 2 


Non-interacting Majorana modes 
(weak coupling to critical 
fluctuations irrelevant) 


Point defect in TSC: Non-trivial 
fixed point in e = 3 — d expansion: 

(xMx(O)) ~ 1/t" with r) w 0.32 


1. Ising magnetic ordering 

of TSC surface: Tricritical Ising point 

2. XY-magnetic/superconducting 
ordering of TI surface: 

central charge c = 2 — > 1 transition. 


d = 3 


Non-interacting Majorana 

modes in TSC and Non-interacting 

Dirac modes/Helical Fermi liquid in 

TI (weak coupling to 

critical fluctuations irrelevant) 


Line-defect in TSC: 

logarithmic corrections to scaling, 

a 'marginal Majorana liquid'. 


1. Ising magnetic ordering 
of TSC surface: 

M = 1 Supersymmetric critical point 

2. Superconducting ordering 

of TI surface (with chemical potential 
tuning): M = 2 Supersymmetric 
critical point. 



Table I: Compilation of results discussed in the main text. The three columns correspond to the three class of problems shown 
schematically in Fig[T] Fig(4]and Fig[5] 




Figure 4: A schematic drawing of a point defect carrying two 
Majorana modes (blue circles) coupled to critical magnetic 
fluctuations at the transition between a non-magnetic two- 
dimensional Dili TSC and a magnetically ordered ordinary 
insulator. 





































7' 



Non-magnetic, TSC 



Magnetic, Trivial SC 



Figure 5: A schematic drawing of a line defect carrying heli- 
cal Majorana modes (blue lines) coupled to critical magnetic 
fluctuations at the transition between a non-magnetic three- 
dimensional Dili TSC and a magnetically ordered ordinary 
insulator. 



modes in certain cases [77|. A general classification of 
defect modes was presented in Ref.jHD]. 

(i) Point defects in two-dimensional TSC at criticality: 



As a first example, consider a point defect in a two- 
dimensional TSC in class Dili that hosts two Majo- 
rana modes \l an d Xi which are coupled to Ising crit- 
ical fluctuations (j> (Figf4|). The total Euclidean ac- 
tion S of the system near the critical point is given by 

S = Sbulk + Sdefect + Sbulk-defect where 



1 



d d xdr( ±(a Al 0) 2 + -0 2 + ^ 4 ()l5) 



Sbulk = 

Sdefect = J dr (xidrXi + xAxz) (16) 

Sbuik-defect = ig J dr Xi(t)x2(t) (/)(t,x± = 0) (17) 

Note that in contrast to the case of surface states cou- 
pled to critical fluctuations studied above, where one 
needs to consider the surface critical exponents for the 
field (j), in this case the relevant exponents are deter- 
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mined by the conventional bulk criticality. This is be- 
cause even though the superconducting order parameter 
vanishes at the location of the topological defect, the de- 
fect does not provide any constraint on the magnetic fluc- 
tuations. This is consistent with the conventional treat- 
ment of Kondo problem impurity and related problems 
[4ll - |46| where bulk criticality determines the fate of im- 
purity spin. The differences with the conventional Kondo 
problem are worth mentioning as well. In particular, the 
'impurity' in our problem is inherently fermionic, and 
satisfies the constraint x\ = x\ = 1 without requiring 
any chemical potential. In contrast, in the Kondo prob- 
lem with, e.g., spin S=l/2 impurity, if one represents the 
spin as S as S = hz^ffz, where 2 is a two-component 
'slave' boson/fermion, one needs to explicitly enforce the 
constraint z^ z = 1 explicitly using a chemical potential 
or by other means |8CH82| . In this sense, our problem has 
features both of a Kondo impurity problem and Yukawa- 
Higgs problem of fermions coupled to bosons [Til]. We 
will perform a perturbative renormalization group (RG) 
in the parameter e = 3 — D, the physically relevant case 
being e = 1. We note that a potential anisotropy term 
/ drcj) 2 (t, x± = 0) will also be generated under RG but 
it is irrelevant since its scaling dimension is unity when 
e = 0. The RG flow for g is (Supplementary Information, 
Sec. B): 



dg 
dl 



2 



(8 - 7r)g 3 
16tt 3 



(18) 



Thus we find a stable infra-red (IR) fixed point at g* = 
-\/87r 3 e/(8 — n). The anomalous dimension of the Majo- 

rana fermion at this fixed point is r\ x — — 2(8-7r) • 
Setting e = 1, one finds (x(t)x(O)) w 1/t 32 . Therefore, 
the Majorana fermion behaves like a zero-dimensional 
strongly coupled system at criticality, similar to the case 
of the Majorana edge states of one-dimensional Dili TSC 
at bulk criticality. A possible physical realization is a dis- 
location in a two dimensional sheet of the organic super- 
conductor (TMTSF) 2 PF 6 as discussed in the following 
section. 

(ii) Line defects in three-dimensional TSC at criticality: 
An analogous problem is that of a line defect embedded 
in a 3 + 1-d TSC in class Dili at the quantum critical- 
ity between the TSC and an Ising magnet (FigJS]) . Away 
from the criticality, the line defects in TSC host free he- 
lical Majorana modes xl an d Xr whose Hamiltonian is 
same as that in Eqn|3J The bulk action is identical to 
Eqn llSI with d — 3 while 



Sbuik-defect = %ig J drdxi Xl(t,Xi)xr(t,Xi) x 

4>{t,xi,xx = 0) 
In addition, a potential term v J dTdxi<j) 2 {r, x\, x± 



0) is also allowed. Unlike the case of point defect embed- 
ded in a two dimensional TSC, a simple power-counting 
shows that in this case, both the coupling g as well as v 
are marginal. The RG flow of g and v is given by (see 
Supplementary Information, sec. C) 



dg 
dl 
dv 
~dl 



31og(2) 3 
16tt 2 9 
3uv 



64tt 2 



Because of the marginality of g, the correlation function 
of the Majorana fermions G(z) = (xr(®)xr( z )) receives 
logarithmic corrections at criticality: 



G(z) 



z(l0g(M)V6 



(19) 



where z = x + it. The correlation function for the 
left-moving Majorana fermion is simple given by the re- 
placement z — > ~z in G. Therefore, the helical Majorana 
metal in the vortex core is a 'marginal Majorana liquid' 
in the quantum critical regime. 

A similar problem can be discussed in the context of 
a dislocation line in a 3dTI, which may carry protected 
one dimensional modes when a weak index is present [77J ■ 
When bulk magnetic ordering sets in via an ising critical 
point, the fate of coupling the ldmodes to these critical 
fluctuations remains to be discussed. Since interactions 
already modify the electronic properties qualitatively, the 
answer depends on further details and will be analyzed 
elsewhere. 



PROSPECTS FOR EXPERIMENTS 

To experimentally realize the predictions above, one 
needs both a symmetry protected topological phase as 
well as proximity to a quantum critical point where a 
symmetry protecting the topology is spontaneously bro- 
ken. Topological insulators like Bi2 Se3 doped with mag- 
netic atoms like Fe, Mn, Cr etc. may undergo a magnetic 
transition either just on the surface or in the bulk, on 
tuning concentration, as recently reported IH 19^ . Cor- 
related systems such as rare earth based half-Heusler ma- 
terials 0, as well as iridium and osmium based oxides 
may also exhibit a topological insulating phase in prox- 
imity to magnetic order 21 , 2^] ■ 

Some of our most interesting results pertain to topolog- 
ical superconductors with time reversal symmetry (DIII- 
TSc). A well known physical realization is the B phase 
of superfluid Hesfi^l- In the solid state context, an in- 
teresting proposal that Cu 2 .Bi2Se3 may be a DIII-TSc 
was recently made(28|. for which there is some intrig uing , 
but as yet not conclusive, experimental indications[29|. 
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Introducing magnetism into this system remains to be 
explored. 

Among the most tunable materials are rare earth inter- 
metallics (heavy fermion systems) and organic materi- 
als, which show a variety of phases and quantum phase 
transitions. A promising set of candidate materials are 
UGe2,URhGe and UCoGe j2(| which show overlap of fer- 
romagnetism and superconductivity in their phase dia- 
grams, and which can be tuned by pressure. If the su- 
perconductivity, which is likely triplet, is topological in 
the paramagnetic phase, then this would be a class of 
candidate materials (37). Finally we mention the quasi- 
one dimensional organic superconductor (TMTSF^PFg 
(30T | is believed to be have equal spin triplet pairing. The 
proposed pairing given the highly one dimensional na- 
ture of the Fermi surface points to a d=l topological 
superconductor in class DHL A two dimensional sheet 
of this material will then feature a weak topological in- 
dex. Interactions with a nearby magnetic phase have 
also been discussed (3lj . These experimental works con- 
tain between them the requisite ingredients to observe 
the physics discussed in this paper, but more work is of 
course required to incorporate all desired properties in a 
single realization. 

CONCLUSIONS AND PERSPECTIVE 

To summarize, we studied the fate of edge modes and 
topological defects in TI/TSC in a wide range of phase 
transitions that do not close the single fermion gap in 
the bulk. Our results are summarized in table I. One 
broad and somewhat surprising lesson we learned is that 
for weak coupling to the critical fluctuations, the edge 
modes that arise due to topological nature of TI/TSC re- 
main well-defined even at the criticality. In some cases, 
such as when the bulk undergoes phase transition in a 
three dimensional TSC, even the qualitative nature of 
the surface states remains unchanged for weak coupling 
to the critical fluctuations. While in other cases, such 
as in the case of Majorana modes at the ends of a one 
dimensional TSC, the localized modes instead show be- 
havior rather different from free Majorana fermions at 
the criticality. We also find evidence for emergent super- 
symmetry at surface phase transitions in TSC/TI in sev- 
eral cases. Measurable consequences of these results for 
experiments as well as numerics were discussed and cur- 
rently available experimental possibilities were reviewed. 
In future it will be interesting to study the effects of dis- 
order on these conclusions, since this would be relevant 
to the case where randomly placed impurities induce a 
magnetic transition. To conclude, the study of critical 
phenomena in topological insulators and superconductors 
lies at the intersection of many conceptual and topical 
themes in condensed matter, such as interacting topo- 
logical phases, bulk and boundary quantum criticality, 



and strongly coupled fermionic liquids, while also open- 
ing more exotic possibilities such as the realization of 
emergent supersymmetry in condensed matter systems. 
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SUPPLEMENTARY INFORMATION 

A: Duality for 1-d TSC coupled to Critical Ising 
Chain for a finite chain 

In the main text, we considered the fate of boundary 
Majorana modes in a semi-infinite one-dimensional time- 
reversal invariant TSC at the phase transition between 
TSC and an Ising ordered state. Here we repete a similar 
analysis for a finite chain that has a pair of Majorana 
modes at both ends. The Hamiltonian for the coupled 
Majorana-Ising system at the bulk criticality is: 

L 

H = -J2 + °f) - w lxiX2 - ioixixa (20) 

i=l 

where we have set the coupling between the Majorana 
modes and the critical chain to unity (the deviation from 
unity turns out to be irrelevant at the stable fixed point, 
via the same argument as for the semi-infinite chain in 
the main text). We also ignore the interaction between 
the Majorana modes at the two ends since it would be 
suppressed exponentially in the system size L. We define 
■1X1X2 = eg and *XiX2 = cr£ +1 : 

L L 

H as written above corresponds to transverse field 
Ising without the transverse field for the first spin 00 
and the last spin ctl+i- Following (59| . let us perform 
the following duality transformation: 

0% = ^ (22) 

of = Mf-iMl V L + 1 > i > 1 (23) 

«H = V L > 1 > (24) 

= Ml+i (25) 

Hamiltonian H in terms of new variables fx is given by 
the standard transverse field Ising model: 

L 
i=0 
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The unequal time correlation function for the Majo- 
rana modes at the two ends are given by (Xa(0)Xa( r )) = 
K(0K(r)) = (/xg(0K(r)) ~ 1/t and ( X ' a (0)x' a (r)) = 
K+i(0)^+i(t)> = K(0)Mi(r)) ~ 1/t for both a = 
1,2. 



modes without any cross-terms. Si, the interaction part 
includes the (j) 4 term as well as Sbuik-edge- Since the <p 4 
part is unaffected by the edge modes, we will only focus 
on the g term as far as Si is concerned. The full partition 
function is given by 



B: Renormalization for point defect in two 
dimensional TSC coupled to critical fluctuations 

Since the bulk is unaffected by the edge modes, we first 
imagine that we have already performed renormalization 
for the bulk degrees of the freedom and are sitting at the 
bulk criticality. The RG equations for Sbuik to 0(e) are 



Z 



D(f><D Xl< D X2< e~ 



D<f>> D X i> D X 2> e~ s ° <3^ 



= z> 



D0 < D X i<D X 2<e- s ° (e- Sl ) 



(33) 



where Zq = J D(j) > D X i > D X2> e s ° and 



du 9u 2 

dT = eu ~ 8^ 

dr 3uA 2 

— = 2r H 

dl 8tt 2 



(e~ Sl )> = ^> J D Xl> D X2> e 



(34) 



3w 
8^ 



(27) 
(28) 



where A is a UV cut-off. This implies that there is 
a stable IR fixed point (with one relevant direction) at 

To perform Wilsonian RG for S e d ge + Sbuik-edge, we 
separate the degrees of freedom into 'slow' and 'fast' i.e. 
we write 



The integration over the fast modes 4>>, X i>, X 2> 
would modify the action for the slow modes 0< , \i< > X2< , 
thus yielding the RG equations so desired. Since we sus- 
pect a fixed point at g* = 0(y/e), a cumulant expansion 
in powers of g is a legitimate option. Therefore, upto 
0(g 3 ), (e _Sl ) may be written as 



/g-SA =e ~((S 1 ) > -^[(sf) > -(S 1 )l]+^[(sf) > -3(sf) > (S 1 ) > +2(S 1 f > }) 

(35) 



x{xq) = 



dk X (k )e lkoxo 

feo|<A/s 

X<(^o) + X>(^o) 



A<|fc |<A/s 



We next write down the contribution from each of the 
dko x(fco)e^^ ee terms in the cumulant expression to the renormal- 
ization of the Sbuik-defect as well as the anomalous di- 
(fl§i)sion of Majorana fermions. 



Similarly, we define, 
(x ,r± = 0) = 



O(g) : (Si). 



dk d d - 1 k± <P(k Q ,k ± )e tkoXo + 

fco|<A/s 

dk d d - 1 k± (f>{k ,k ± )e lk ° X0 

A<\k \<A/s 

= Mr)+0>(r) (30) 

where r = (xo,rj_). Note, in particular, that the per- 
pendicular momenta k^ are being integrated over the 
whole range in both <j> < and </>>. This is legitimate be- 
cause ~ do not generate any singularities for the 
edge action as long as the modes near fco ~ are never 
integrated over, which is indeed true owing to the above 
separation. 

The action S in terms of new variables <^> , <p < , X> , X< 

is 



(Si)> = i 9 / dx Q Xi<(^o)X2<(a;o)0<(a:o,^± = 0) (36) 



Consider the elementary coarse-graining transforma- 
tion x' = Xq/s. Under this transformation, Xi<i x o) = 
txXi(%'o,), X2<(x ) = £, xX 2(x' ,) and 0<(x o ,f± = 0) = 
^(^(xq, x± = 0). £ x and ^ are related to the scaling 
dimensions A x , A^ of the fields £ and X for d = 4 — e 



via 



= s -Vx/2 an d ^ = S -A* = s -(d-i+^)/2_ 



Re-expressing (Si) > in terms of rescaled variables, 



(Si) = igs e/2 I dx' dx' x X T (x' , x'{)^sX(x'q, x'i) <t>i^o, x 'n %± = { 



S — S^~ + Sq + Si 



(31) 



where we have dropped the contribution from r/^ since 
it is of order 0(e 2 ) (we will soon see that r\ x = 0(e)). 
Therefore, at this order, g(s) = gs e ^ 2 ~ Vx where e = 3 — d. 



where So denotes the Gaussian part of the full action 
and which, therefore, separates into the slow and fast 



0(g 2 



l r 



<Si>> - (Si) 
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These terms renormalize the propagator of the edge 
X fermions leading to an anomalous dimension for the 
fermions as shown in diagram. Denoting this contribu- 
tion by I x , 



d ko d d fej 



C: Renormalization for line defect in three 
dimensional TSC coupled to critical fluctuations 

The total Euclidean action of the system near the crit- 
ical point is given by: 



2 ( 27 4 7 <|fc_ L |<A 7A/ s <|fc„|<A Oo+R))0o+ &i+y 



S — Sbulk + Sdefect + Sbulk-defect 



(47) 



9 



16tt 2 



-po log(s) 
Hence rj x = 



^ggwhere Sbak is same as in the case of point defect with 
5 = 3 and 



S,i, 



efect 



0(g 3 ): -[(S 3 1 ) > -3(S 2 ) > (S 1 ) > +2(S 1 ) 



Only the term (S?) > contributes to the renormaliza- 
tion of the interaction g at this order as depicted in FiglHl 
Let us call these contribution AS g . 



AS g = x(ig) 3 / dx Xi<(x )X2<(x ) (f><{x , x± = 0) X (240) 



where 



J dx dxi (xrOtXr + XlOtXl + iXRdxXR - *XiAXi) 
dx dx x x t 7m^X (48) 



where the matrices 70 = t x , 71 = r y act on the vector 

X = [xrXl] T and x T = X T 7o- Tne coupling S bu i k ^ de f ect 
between the critical bulk and boundary CFT is given by 

Sbulk— defect 

2ig / dx dxi Xl(xo,xi)xr(xq,xi) <j>(xo,xi,x± = 0) 



1 



(2tt)* 

1 



dpo 



A/s<|p|<A 



d 3 p 



o<| P J<A Po(p'o+Pl + r 



ig / dx dxi x T (x , xihsX^o, x{) (j>(x , xi, x± =6) 



2^ X (1 - n/4) X V 



Using this, one finds 



t \ e/2— J 
l{S) = S ' 16 



where 75 = r z = —17071 ■ In addition, one generates 
^ ' a potential anisotropy term v J dxodxi<p 2 (xo, x±, x± = 0) 
localized near the defect. 

Again the bulk criticality is unaffected by the bound- 
ary. Similar to the case of point defect, we write 



1 - (s - 1 



9 3 ) (43) 



X(r) 



Writing s = e , the RG flow for g becomes 

dg = eg_ (8 - n)g 3 
dl 2 16tt 3 



d 2 k X {k)e lk - F + 

\k\<A/s 

X<{r) + X>(r) 



d 2 k x{k)e 



ik.f 



A<|fc|<A/s 



(49) 



(44) 



Thus, we find a stable IR fixed point at g* = 
^87r 3 e/(8- tt). 

The anomalous dimension of the Majorana fermion is 
given by 



where r= (xq,Xi) and 
fi(r, r±=0) = 



d 2 kd d ~ 2 k ± <j>(k,k ± )e ik - p 

k\<A/s 

d 2 kd d - 2 k ± <f>{k,k ± )e ik - p 

A<|fc|<A/s 

4><(f) + <P>(r) 



(50) 



9 



*2 



The full partition function is given by 



16tt 2 
en 



2(8 -it) 

Setting e = 1, one finds r\ x = 2 (%--n) 



(45) 

(46) Z 



0.32. 



= J D(/) < Dx<Dx< e ~ S ° J D^DxyDxle" 3 " <5f) 
= Z> J D^Dx^xle- 3 " (e- s >) > (52) 
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where Zq = J DfyyDxyDx^e s » and 




<e- Sl )> =-g>J D^DxyDxle-^e-^ 
Upto 0(g 3 ), (e~ Sl ) > may be written as 

renormalization ot vertex g at U(g") tor trie problem ot a 
/ e -Si\ _ e -(<Si> > -^[(s?) > -<Si)?,J+i[(Si ! ) > -3(5i ! ) > <Si> > +2< mll g > |j^f ect COU pi ec i to critical fluctuations. The wavy line 



Figure 6: The Feynman diagram that contributes to the 
renormalization of vertex g at 0(g 3 ) for the problem of a 



(53) 

We next write down the contributions from each of the 
three terms in the cumulant expression. 



represents the propagator H^(p, p±) = +r for the bulk 

<j> field while the straight line corresponds to the propagator 



lb 



for the edge fermions. 



0(g) ■■ (Si} : 



(Si) > = ig / dxodxt X T <{xo,xi)~f S X<(x , %i) <t><{xo, Xi, x± = 0) 



dv 



64tt 2 



(55) 



Consider the elementary coarse-graining transforma- 
tion Xq = xq/s, x'i = xi/s. Under this transforma- 
tion, x<(xo,xi) = CxXC^oj^'i) and <f><(xo,Xi,x± = 0) = 
^^{Xq j X^ , X_\_ = 0). £ x and are related to the scaling 
dimensions A x , of the fields £ and x a t the free fixed 



0(g 3 ): -[(S 3 ) > -3(S 2 ) > (S 1 ) > +2{S 1 ) 



point via £v 



-A 



x = s -(i/2+ nx /2) and £ = s - 



A, 



s ^ d 1+, '<*)/ 2 . Re-expressing (Si) > in terms of rescaled 
variables. 



(Sf') > contributes to the renormalization of the inter- 
action g at this order as depicted in FiglBJ 



AS g = (ig) 3 / dx dxi x T <(xq,xi)^sX<(xo,xi)4' < (xq,xi 1 x± 



(Si 



igs 



£/2-Vx 



dx' dx[ x T {x' , x'lhsxix'o, x[) </>(x' Q , x[, fwfaeift) 



Therefore, at this order, g(s) = gs t l 2 Vx where e = 
3-d. 



d 2 p 



0{g\uv): -l[( S f) > -(S 1 ) 



These terms renormalize the propagator of the edge 
X fermions leading to an anomalous dimension for the 
fermions as shown in diagram. Denoting this contribu- 
tion by I x 



A/s<|p|<A 

log(2) a-1 
~ 8tt 2 s 

Collecting everything, 



P± 



0<\p±\<A 



p 2 (p 2 +p]_ + r*) 



(56) 



e/2 2_ 

S 16n 



log(2)s-l 3 



d 2 k\\ d 2 k±_ 



(2tt) 4 y <|fc_ L |<A A/ s <|fc|,|<A 0>7m + Pm7m)(^ 2 + k ±) 



Writing s = e dl , the RG flow for g becomes 



167T 



;Pn1n log(s) log(2) 



(54) 



dg_ = t9__ 31og(2 ) 
dZ 2 16tt 2 



(57) 



(58) 



Hence % = T |^. 

One also generates renormalization of v: 

AS = <f) 2 < (x ,x 1 ,x± = 0) x 



For the problem of line defect in a three dimensional 
TSC, e = 0, hence 



1 



dfc|| d d k± 

( 2?r ) 4 Jo<|fci|<A Jp<|fe|||<A ( fc2 + fc l) 2 

This leads to 



dg _ 31og(2) 3 

2 ^ 



16tt 2 



(59) 



To the leading logarithmic correction, the renormaliza- 
tion group equation for the correlation function G[z, g) = 
(Xr($)Xr{z)) is gi ven b y 
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Figure 7: The Feynman diagram that contributes to the self- 
energy of Majorana fermions for the problem of a line defect 
coupled to critical fluctuations. 



rlC r)C 
dlog(z) dg 

where P{g) = = -Ag* with B = ^ and 

A = 3 ■ Solving the above differential equation, one 
finds 



G(z) = - } B (61) 

Z(l0g(|2|)2* 

= z(log(|z|)i/6 ^ 
For the left-moving fermions, one simply replaces z — > 
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